Draft version February 5, 2008 

Preprint typeset using LSTgX style emulateapj v. 10/09/06 



CONSTRAINTS ON THERMAL X-RAY RADIATION FROM SAX J1808.4-3658 AND IMPLICATIONS FOR NEUTRON 

STAR NEUTRINO EMISSION 1 

C. O. Heinke 2 - 3 , P. G. Jonker 45 - 6 , R. Wunands 7 , and R. E. Taam 2 

Draft version February 5, 2008 

ABSTRACT 

Thermal X-ray radiation from neutron star soft X-ray transients in quiescence provides the strongest con- 
straints on the cooling rates of neutron stars, and thus on the interior composition and properties of matter 
in the cores of neutron stars. We analyze new (2006) and archival (2001) XMM-Newton observations of the 
accreting millisecond pulsar SAX J1808. 4-3658 in quiescence, which provide the most stringent constraints 
to date. The X-ray spectrum of SAX J1808. 4-3658 in the 2006 observation is consistent with a power-law of 
photon index 1.83 ±0.17, without requiring the presence of a blackbody-like component from a neutron star 
atmosphere. Our 2006 observation shows a slightly lower 0.5-10 keV X-ray luminosity, at a level of 68^)3% 
that inferred from the 2001 observation. Simultaneous fitting of all available XMM data allows a constraint on 
the quiescent neutron star (0.01-10 keV) luminosity of L^s < 1.1 x 10 31 erg s" 1 . This limit excludes some cur- 
rent models of neutrino emission mediated by pion condensates, and provides further evidence for additional 
cooling processes, such as neutrino emission via direct Urea processes involving nucleons and/or hyperons, in 
the cores of massive neutron stars. 

Subject headings: binaries : X-rays — dense matter — neutrinos — stars: neutron 



1. INTRODUCTION 

The X-ray transient SAX J1808.4-3658 (hereafter 1808) 
has provided many fundamental breakthroughs in the study 
of accreting neutron stars (NSs). It was discovered in 1996 
by BeppoSAX's Wide Field Cameras, and type I X-ray bursts 
were seen, ide ntifying it as an accret i ng NS and constraining 
the d istance din 'tZand et alJ 1 1998t iGalloway & Cum ming 
2006). Coherent millisecond X-ray pulsations, the first 
discovered in accreting systems, were identified during an 
outburst using the Rossi X-r ay Timing Explorer (RXTE; 
Wiina nds & van der K lis 1998). Burst oscillations have also 
been seen at 1808's 401 Hz spin frequency, confirming that 
thermonuclear burst oscillations in low-mass X-ray binaries 
(LMXBs) represent the spin pe riod of the NS (lin't Zand et al.1 
2001; Chatoa b^t7et^l2003h . A pair of kilohertz quasiperi- 
odic oscillations (QPOs) were seen from 1808, with a fre- 
quency difference equal to one-half of the spin period, 
forcing a revision of the most popular models for QPOs 
(Wiin ands et al.l l2003h . Optical observations of the brown 
dwarf companion while 1808 was in quiescence showed a si- 
nusoida l optical modulati on attributed to heating of the com- 
panion dHomer et alj200ll) . However, the required irradiating 
luminosity is larger than the available X-ray luminosity, giv- 
ing rise to speculatio n that a radio pulsar mechanism is activ e 
during quiescence dBurderi et alJ2003llCampana et alf 2004). 

In addition to these discoveries, 1808 has provided one 
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of the lowest quiescent thermal luminosities yet measured 
from any accreting NS dCampana et al.ll2.002i: IWiinands et all 
12002b . along with 1H1905+000 dJonker et alj |2006h. Tran^ 
siently accreting NSs in quiescence are usually seen to have 
soft, blackbody-like X-ray spectra, often accompanied by 
a harder X-ray co mponent generally fit by a power-law of 
photon index 1-2 dCampana et al.lfl998b . The harder com- 
ponent is of unknown origin; an effect of continued accre- 
tion, or a shock from a pulsar wind have been suggested 
(Cam pana et al.l [T998b . The blackbody-like component is 
generally understood as the radiation of heat from the NS sur- 
face. This heat is produced by deep crustal heating during 
accretion, and is radiated by the crust on a timescale of 10 4 
years, producing a steady quiescent thermal NS luminosity 
Brown et alj|1998t ICampana et all [19981: lHaensel & Zdunikl 
11990b . Measurement of the blackbody-like component is par- 
ticularly important because it indirectly constrains the interior 
structure of NSs. 

What fraction of this heat escapes the NS as neutrinos rather 
than photons depends on the physical conditions (i.e. com- 
position, density, and pressure) of the NS interior. If the 
outburst history (fluence, recurrence time) and distance of 
a NS are reasonably well-known, then the determination of 
the quiescent thermal NS luminosity constrains the neutrino 
vs. photon emission, and thus models for the NS in terior 
dYakovlev & Pethickll2004t iLevenfish & Haenselll2006h . For 
example, the transient Cen X-4 has been identified as hav- 
ing a rather low quiescent X-ray luminosity compared to deep 
crustal heating predictions. This suggests that Cen X-4 has 
enhanced neutrino emission, prod uced in the high density 
core of a relatively high mass NS dColpietalJl200lh . Many 
other LMXBs have also shown low quiescent thermal X-ray 
luminosities, indicating either enhance d neutrino emission or 
extremely long quie scent intervals (e.g.lWiinands et aT1l2001b 
Uonker et al]|2004bt rTomsick et alj|2004t Uonker et alj|2006h . 
The coolest of these provide the strongest constraints to date 
on neutrino cooling from NS cores, as a broader range of cool- 
ing rates is necessary to explain the data than for young cool- 
ing pulsars dPage et al.ll2004l;lYakovlev & Pethickll2004}) . 
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1808 is a particularly interesting system due to its 
known distance and constr aine d mass transf er rate 
dBildsten & Chakrabartvl 1200 lb . ICampana et ail d2002b 
observed 1808 in quiescence with XMM in 2001, finding a 
low luminosity (L z (0.5-10keV)= 5 x 10 31 ergs s" 1 ,ford = 2.5 
kpc) and a relatively hard spectrum, with less than 10% of the 
X-ray flux attributable to a possible blackbody-like compo- 
nent. We have obtained a deeper XMM observation of 1808 
in quiescence in 2006 to place more stringent constraints on 
neutron star cooling processes. 

2. DATA REDUCTION 

We observe d 1808 on September 14, 2006 
Jobs-id 0400230 101') for 54 ksec wi th XMM's EPIC 
camera, using two MOS CCD detectors (Turner et al.ll2001h 
with medium filters and one pn CCD detector ( Striid er et al.l 
2001) with a thin filter. We also downloaded the March 
24, 2001 XMM observa tion (obs-id 0064940101, reported 
by ICampana etaT1l2002h from the HEASARC archive. All 
data were reduced using FTOOLS and SAS version 7.0.0. 
We used only the MOS data from 2001, since the pn data 
were taken in timing mode, and the target was too faint to be 
detected in this mode. Intervals of flaring background were 
excluded by excluding times when the total MOS count rate 
exceeded 5 0.2-12 keV counts per second, and times when 
the total pn count rate exceeded 50 0.2-12 keV counts per 
second. This left 36.6, 47.0, 48.1, and 39.3 ksec in the 2001 
MOS data, the 2006 MOS1 data, the 2006 MOS2 data, and 
the 2006 pn data, respectively. Event grades higher than 12 
were also excluded. We extracted spectra from a 10" circle 
around the position of 1808, correcting the fluxes for the 
fraction of photons collected within the extraction radius, 
and combining the pairs of simultaneous MOS spectra and 
responses using FTOOLS. We generated response and effec- 
tive area files using the SAS tasks rmfgen and arfgen, and 
produced background spectra from 90 " circular source-free 
regions on the same CCD. The spectra were grouped to 
>15 counts per bin for the MOS data, and >30 counts per 
bin for the pn data (other choices gave similar results). To 
assess variability within the 2006 and 2001 observations, 
background-subtracted lightcurves were produced within 
SAS and analyzed using HEASARC's XRONOS software. 
KS and \ 2 tests on the l ast 15 ksec of 0.2-12 keV pn data 
(unaffected by background flaring) revealed no evidence of 
variability. 

3. SPECTRAL ANALYSIS 

Our fitting includes photoelectric absorption (XSPEC 
model phabs), with a hydrogen column density, Ng, fixed 
at the value derived from observations in outburst (1.3 x 
10 21 cm" 2 ; note this is equal to the value derived from 
iDickey & Lockmanl [l990). We also tested models with pho- 
toelectric absorption as a free parameter, in all cases finding 
Nh consistent with the outburst value. Quoted errors are at 
90% confidence. 

We fit the combined 2001 MOS spectrum to a power- 
law model, finding T = 1.72 ±0.28 (see Table 1). Fits 
using on ly a hydrogen-atmos phere model, the NSATMOS 
model of iHeinke et al] (l2006ai) (similar to the NSA model of 
IZavlinetal.ll 1996b gave poor fits (x 2 > 4.8). We then per- 
formed fits with NSATMOS plus a power-law, fixing the true 
NS radius to 10 km, the gravitational mass to 1 .4 M , and the 
distance to 3.5 kpc (Gallowav & Cumming 2006). No thermal 
component is required, but a thermal component with kT <42 
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FIG. 1. — Top panel: XMM X-ray spectra (data and best-fit power-law 
model) of SAX J 1808.4-3658. Solid line with crosses: 2006 pn data, model. 
Dashed line with triangles: 2006 MOS data, model. Solid line with circles: 
2001 MOS data, model. Bottom panel: residuals to the fit in units of \ 2 . 

eV (90% confidence 8 ) is permitted, thus placing a limit on the 
NS's thermal 0.01-10 keV (essentially bolometric) luminos- 
ity of Lns < 2.4 x 1 31 erg s" 1 . The inclination of t his system 
is known to be low (Bildsten & Chakrabarty 2001), and it is 
rare for L MXBs to show higher Nh in quiescence than out- 
burst (e.g. Jonker & Nelemans 2004]), so we regard a higher 
Nh as very unlikely. The total 0.5-10 keV unabsorbed lumi- 
nosity is L x = 7.6+J g x 10 31 erg s" 1 . 

For the 2006 data, we find a similar spectral shape, and 
therefore fit similar models to the pn and MOS data simul- 
taneously. A simple power-law fits the data adequately, with 
a photon index of T = 1.83 ±0.17 and an unabsorbed Lx (0.5- 
10keV)=5.2±0.7 x 10 31 ergs" 1 . The 2006 flux appears less 
than the 2001 flux. We test this by fitting the spectra simul- 
taneously and tying their power-law slopes together, finding 
that the 2001 0.5-10 keV unabsorbed flux is higher at 97% 
confidence; the 2001 flux is 1.28^ (90% conf.) that of the 
2006 observation. If the power-law slopes are allowed to vary 
(1.61 ±0.3 for 2001, 1.83 ±0.17 for 2006), then the best-fit 
flux ratio is 1.47^. 

Fitting the 2006 and 2001 data allows a tighter constraint 
on the presence of a NS atmosphere component than the 2001 
data alone, requiring a NS kT < 34 eV and a thermal 0.01- 
10 keV NS luminosity L NS < 1.1 x 10 31 erg s" 1 . Choosing 
a NS radius of 12 km, or a mass of 2.0 Mq, varies this 
constraint by only 3% . The rather tight distance limits of 
lGalloway~& Cumming (2006) (3.5 ±0.1 kpc) produce only 
a 6% uncertainty. Allowing the Nh to float freely permits a 
thermal 0.01-10 keV NS luminosity L NS < 1.0 x 10 32 ergs 
s _1 (for N H = 1.7 X 10 21 cm" 2 ). 

4. RAMIFICATIONS 

We have estimated the time-averaged mass transfer rates for 
1808, and several other transient LMXBs (Aql X-l, Cen X- 
4, 4U 1608-52, KS 173 1-260, RX 1709-2639, MXB 1659- 
29, XTE 2123-058, SAX 1810.8-2609, and those in Terzan 
5 and NGC 6440), from the RXTE All-Sky Monitor (ASM) 
record (1996 to Nov. 2006) under the assumption that the 
time-averaged mass accretion rate over the last 10 years re- 

8 kT is the redshifted temperature, or kToo 
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FIG. 2. — Cooling curves for various NS interior neutrino emission 
scenarios, compared with measurements (or 95% conf. upper limits) of 
the quiescent 0.01-10 keV NS luminosity and time-averaged mass trans- 
fer rate for several N S transients (see Table 2). Cooling curves taken from 
lYakovle v & Pethick 12004); the dotted curve represents a low-mass NS, 
while the lower curves represent high-mass NSs with kaon or pion conden- 
sates, or direct Urea (Durca) processes mediated by nucleons or hyperons. 
Limits on the quiescent NS luminosity of SAX J1808. 4-3658 are given for 
the 2001 and 2006 observations. The effect of a distance error as large as a 
factor 1.5 is also indicated (upper left). 

fleets the time-averaged mass transfer rate (Table 2). We use 
PIMMS and a power-law of photon index 2 to convert the 
ASM countrates during outbursts into 0.1-20 keV fluxes 9 . 
This is, of course, a rough approximation, as the spectral 
shapes of LMXBs in outburst vary substantially. Additional 
sources of potential error include poor ASM time coverage 
of some outbursts, uncertainty in the NS mass and radius (af- 
fecting the energy released per accreted gram, and thus the 
conversion from Lx to mass accretion rate), variability in the 
mass transfer rate, and uncertain distances (which will equally 
affect the quiescent luminosity). We plot an arbitrary uncer- 
tainty of 50% in both mass transfer rate and quiescent lu- 
minosity for each point in Fig. 2. For Cen X-4 we use the 
lowest measured quiescent luminosity, and the mass transfer 
rate limit inferred if Cen X-4 undergoes outburs ts every 40 
years with a fiuence similar to its 1969 outburst (Che n et alj 
119971) . The NS component flux for Aquila X-l is some- 
what uncertain and pos sibly variable (Rutledge et al.l 120021 : 
ICampana & Steiial l2003). We assume that all outbursts from 
NGC 6440 since 1971 have been detected. For KS 1731-260 
we assume that the average flux seen with RXTE-ASM dur- 
ing outburst was the average flux during the entire 12.5 year 
outburst. For KS 1731-260 and the transient in Terzan 1 (for 
which we take a 12-year outburst) we take a minimum recur- 
rence time of 30 years. 

For 1808 we derive a time-averaged mass transfer rate of 
1.0 x 10 -11 M /year, an excellent match to the prediction 
of general relativity of 0.95 x 10- n (M 2 /0.05 AT©) M Q /year 

9 We have verified that this conversion is correct to within 50% for out- 
bursts of the transients EXO 1745-245 and Aquila X-l. 



(Bilds ten & Chakrabartvll2001h . We note that the true mass 
transfer rate cannot be less than 7 x 10~ 12 M /year for a NS 
mass >1.4 M under the assumption of an index n = -1/3 
for the donor's mass-radius relation, or less than 3.5 x 10~ 12 
M /year for an index n = 1 mass-radius relation. It is un- 
clear whether the entropy of the donor can be ma intained by 
the low quiescent luminos ity of the neutron star dHomer et alj 
1200 U iBurderi et ail [2003): we plan future observations and 
modeling to address this issue. 

We have plotted the cooling curves calculated by 
lYakovlev & Pethickl (120041) for a variety of models in Fig. 2. 
Low-mass NSs will cool slowl y (do tted line in Fig. 2); in the 
model of Yako vlev & Pet hick (2004), only through photon 
emission and neutron-neutron neutrino bremsstrahlung pro- 
cesses, while modified Urea neutrino emission is suppressed 
by proton superfluidity. Other slow cooling models (invok- 
ing, e.g., neutrino emis sion through Co oper pair formation) 
give similar results (e.g. Page et al. 2004). Higher-mass NSs 
should have higher central densities, sufficient to promote 
more rapid direct Urea neutrino cooling processes involv- 
ing nucleons and/or hyperons, or direct Urca-like processes 
mediated by pions, kaons, or quark matter, in their cores. 
Medium-mass models can produce intermediate cooling rates 
if proton superconductivity is important at low densities, as 
its decay at moderate densities can allo w a smooth transition 
between fast and slow co oling rates (Yak ovlev et al.l 120031 : 
iLevenfish & Haensell2006l) . Thus NSs of different masses 
should lie between the top curve and one of the lower curves 
in Fig. 2, where the lower curve is the maximum neutrino 
cooling curve. We note a possible trend that NSs with low 
mass transfer rates seem to have particularly low quiescent 
luminosities, well below the "standard cooling" predictions. 
This might be explained through binary evolution; NSs with 
low mass transfer rates may be very old systems, which may 
have accreted significant mass. In enhanced neutrino emis- 
sion scenarios, these massive NSs would then have higher 
neutrino and lower photon luminosities than younger systems. 

Our constraint on the quiescent thermal 0.01-10 keV NS 
luminosity of 1808 from the 2006 observations thus seems 
to rule out some models of direct Urea neutrino emission 
via pion condensates, favoring direct Urea processes involv- 
ing nucleons and/or hyperons. An e xtremely large distance 
uncer tainty of 50% (see Fig. 2; cf. Gall oway & Cummingl 
2006) would be required to bring 1808's thermal luminos- 
ity up to the pion condensate predictions. The 2001 obser- 
vation, by itself, rules out some models of direct Urea neu- 
trino emission from kaon condensates. Other modelers of NS 
coolin g have suggested that me dium effects (|Blaschke et alj 
120041) or diquark condensates (Grigor ian et al.l 120051) could 
provide a wide range of NS cooling rates. These models 
may also be sufficient to explain the data on 1808 presented 
here. Our results agree with the principal conclusions of , e.g., 
Yakovlev & Pethickl d2004l) . lLevenfish & Haensell (120061) . and 
Page et alJd2006l) . and provide a firmer observational basis for 
future studies. 
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TABLE 1 

Spectral Fits to SAX J1808.4-3658 



Epoch 


N H 
(10 22 cnr 2 ) 


r 


X 2 /dof 


Lx 
(erg s" 1 ) 


kT 
(eV) 


Lms 
(erg s~') 


2001 

2006 
2001 & 2006 
2001 & 2006 


(0.13) 
(0.13) 
(0.13) 
0.15±0.04 


1.61 ±0.3 
1.83±0.17 
1.83±0.16 


0.51/9 
0.86/45 
0.79/55 
0.78/54 


7.6t\ J 5 x 10 31 
5.2±0.7 x 10 31 
5.2±0.7 x 10 31 
5.2±1.0x 10 31 


<42 

< 35 

< 34 
<61 


<2.4x 10 31 

< 1.2 X 10 31 

< 1.1 x 10 31 

< 1.0 x 10 32 



NOTE. — Spectral fits with power-law plus NSATMOS model to SAX J1808.4-3658. Errors are 90% confidence for a single parameter. Nh is 
held fixed in the first three rows. Lx for 0.5-10 keV range, L N s for 0.01-10 keV. 



TABLE 2 

Luminosities and Mass Transfer Rates 



Source N H kT D Outbursts Years M L NS Refs 

(10 22 cm" 2 ) (eV) (kpc) (M Q yr 1 ) (ergs" 1 ) 



AqlX-1 


4.2 x 10 21 


~94 


5 


8 


10.7 


4 X 10-'° 


5.3 x 10 33 


1,2,3,4 


Cen X-4 


5.5 x 10 20 


76 


1.2 






<3.3x 10-" 


4.8 x 10 32 


5,3 


4U 1608-522 


8 x 10 21 


170 


3.6 


4 


10.7 


3.6 x i(r 10 


5.3 x 10 33 


6,3,4 


KS 1731-260 


1.3 x 10 22 


70 


7 


1 


30 


< 1.5 x 10- 9 


5 x 10 32 


7.4 


MXB 1659-29 


2.0 x 10 21 


55 


-10? 


2 


10.7 


1.7 x 10- 10 


2.0 x 10 32 


7.4 


EXO 1747-214 


4 x 10 21 


< 63 


< 11 






< 3 x 10-" 


< 7 x 10 31 


8 


Terzan 5 


1.2 x 10 22 


< 131 


8.7 


2 


10.7 


3 x 10- 10 


< 2.1 x 10 33 


9,10,4 


NGC 6440 


7 x 10 21 


87 


8.5 


3 


35 


1.8 x 10" 10 


3.4 x 10 32 


11,4 


Terzan 1 


1.4 x 10 22 


74 


5.2 






< 1.5 x lO" 10 


< 1.1 x 10 33 


12 


XTE 2123-058 


6 x 10 20 


<66 


8.5 


1 


10.7 


<2.3 x 10-" 


< 1.4 x 10 32 


3,4 


SAXJ1810.8-2609 


3.3 x 10 21 


< 72 


4.9 


1 


10.7 


< 1.5 x 10-" 


< 2.0 x 10 32 


13,3,4 


RX J1709-2639 


4.4 x 10 21 


122 


8.8 


2 


10.7 


1.8 x ur 10 


2.2 x 10 33 


14,15,4 


1H 1905+000 


1.9 x 10 21 


<50 


10 






< 1.1 x 10-'° 


<4.8x 10 31 


16,15 


SAXJ1808.4-3658 


1.3 x 10 21 


< 34 


3.5 


5 


10.7 


1.0 x 10" 11 


< 1.1 x 10 31 


17,4,15 



NOTE. — Estimates of quiescent thermal luminosities from neutron star transients, and mass transfer rates (inferred from RXTE ASM observations for systems with 
RXTE-era outbursts). Quiescent thermal luminosities are computed for the unabsorbed NS component in the 0.01-10 keV range. Outbursts and years columns give the 
n umber of outbursts and the time baseline used to comp ut e M, if this calculation was performed in this work (indicated by referring t o reference 4). Reference s as follows: 
LIRutledge et alj<2001al) 2:ICampana & StellH <2003h. 3 : ITomsick et~aT] <2004T). 4: M ass transfer rate computed in t his work, 5: Rutledge et al. 1 2001b), 6: Rutledg e~eTaT] 
11999I). 7: ICackett et all <2006bl). 8: ITomsick et all <2005l) 79i Twiinands et all 420051) . 10: IHeinke et al] <2006bl) ll:ICack ett et al. 1 2005), 12: Cackett et al. 1 2006j) ,13: 
Honker et all <2004bl) . 14: Honker et all )2004al) . 15: Quiescent bolometric luminosity computed in this work, 16: IJonker et all {2006b . 17: lGallowav & Cumm ing 1 2006). 



